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Abstract 

We analyse some acoustic emission time series obtained from a lathe machining process. Considering the 
dynamic evolution of the process we apply two classes of well known stationary stochastic time series models. 
We apply a preliminary root mean square (RMS) transformation followed by an ARMA analysis; results thereof 
are mainly related to the description of the continuous part (plastic deformation) of the signal. An analysis of 
acoustic emission, as some previous works show, may also be performed with the scope of understanding the 
evolution of the ageing process that causes the degradation of the working tools. Once the importance of the 
discrete part of the acoustic emission signals (i.e. isolated amplitude bursts) in the ageing process is understood, 
we apply a stochastic analysis based on point processes waiting times between bursts and to identify a parameter 
with which to characterise the wear level of the working tool. A WeibuU distribution seems to adequately describe 
the waiting times distribution. 

1 Introduction 
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When a solid body is subjected to a varying stress, acoustic waves (i.e. pressure waves propagating inside of the 
body), often reaching very high frequencies, are generated. Examples of this phenomenon are the noise produced 
by hitting a metal block with a hammer or the creaking of a wooden floor; the waves propagate inside and on the 
surface of the materials before dissipating in the surrounding gaseous medium. The process of generating acoustic 
waves in stressed materials is called acoustic emission (AE). The acoustic emission can be recorded by means of 
a transducer (i.e. a sensor in contact with the solid body which transforms the elastic waves into electric signals). 
Usually the emission of acoustic waves may also be associated with microfractures inside the solid or, in general, a 
degradation of its condition. Therefore an analysis of the acoustic emissions can reveal the level of degradation of 
the solid. 

In particular we note that the study of acoustic emission is developing in the field of tool condition monitoring 
(TCM) where it is important, for example to avoid damaging machinery and to maximise productive capacity. 
AE analysis is an easily implemented and economical technique and it allows real-time monitoring of working tool 
conditions. 

Oth er statistical analyses have been previously cond ucted using the same experimental setup described in this 
article ( Farrellv et al.l . l2004t IPontuale et al.l l2003i 12001 ). Results demonstrated that acoustic emission analysis is 
particularly relevant for the study of the ageing of work tools and the authors elucidated methods for the elicitation 
and reconstruction of the pdf of the root mean square amplitude signal. 

In this study we analyse some sp ecific high frequency acoustic em issions by means of time series statistical- 
probabilistic models: ARM A models (|Box et al.l , l2004t iBattaglial 120071 ) and models for point processes (JBattaglia , 
20071 ; I Cox and Millerl . Il965l) . An objective of this study is to assess their suitability for explaining AE phenomena, 
and to estabilish if it might be possible to implement a tool wear monitoring algorithm. 

The experimental setup (presented in figure [T]) consisted of a mechanical lathe working on stainless steel bars 
and a transducer with which the resulting acoustic emissions were recorded. This signal was preamplified and 
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Figure 1: The acoustic emissions were recorded by means of a transducer, preamplified, filtered, digitized and then 
stored. 



filtered by a band-pass filter to isolate the frequencies of interest. The digitization was performed by means of a 
digital oscilloscope with a sample frequency of /o = 2.5 Mhz. The cutting speed ranged from 0.5 to 1 ms^^ and 
the cutting depth was 2 mm. 

Three types of working tools were used for the analyses: new, partially and totally worn. For the new and 
totally worn tools we conducted one acquisition with each of eight different tools. For the partially worn tools only 
four acquisistions were conducted. For each acquisition we recorded 15 AE time series, each composed of 40960 
consecutive points (i.e. 614400 points for each acquisition). The time duration of a sin gle time series is 16.38 ms. 

Al l analyses were conducted with the R environment for statistical computing (JR Development Core Team 
(20081), http://www.R-project.orgl, a very powerful programming environment mainly used for statistical data 



analysis and modelling. 



2 Data collected and preliminary analyses 

When we take a close look at the collected time series that result from the experiments, it is possible to distinguish 
two different and superimposed parts of the signal. As we can see in the lower part of figure [21 there is a continuous 
part characterized by a relatively constant variability (essentially due to plastic deformation of the material) and a 
more interesting part composed of bursts of different and high amplitudes (usually associated with microfractures 
in the tool or with material splinters striking it). 

The natural time scale of the cutting process and the short duration of e ach time serie s suggest that each time 
series will be effectively stationary; this is verified by the Dickey-Fuller test ([Greenel . l2003i par. 18.3.3). 

Because of the large amount of data (due to the high sample rate necessary to capture the highest fre- 
quencies), it is necessary to operate a data reduction by means of a transform that has some physical meaning 
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Figure 2: A typical AE signal with its RMS transform above with a threshold at level 20. 



(JKannatev-Asibu and Dornfeldl . Il982l ) 



Definition 2.1 (RMS Values) Let Xg, s ^ N be a time series taking values in 
in some interval over which the RMS yt is to be calculated. We define 
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If a time series is of length N (multiple of T), then the RMS series j/t G M has N /T points, each proportional 
to the acoustic energy emitted in the interval {tT — T, tT\ . Our choice of T is motivated by examination of the AE 
spectral density function a;(///o), which is naturally separable into lower and higher frequency regions by an almost 
zero density region centred at ///o ^ 5 x 10~'^. Furthermore, analysis of x{f / fg) on a moving window demonstrates 
that the low frequency region is largely constant while the higher part depends critically on the presence or absence 
of bursts in the window. 

As bursts are typically separated by the 0(10*) samples, we conclude that the lower part of the spectrum is due 
to the process generating the series of bursts, whereas the higher part is due both to the dynamics of individual 
bursts and to the structure of the continuous part of the AE signal. In order to obtain a meaningful RMS, that is, 
to reduce the data to a small number of points characterising the AE noise level, we chose T = 100 corresponding 
to the spectral gap at ///o = 5 x 10"'^ (equal to a physical interval of 40^s). 



3 Results for ARMA models 

We now subject this RMS series yt to analysis to determine if some feature of the resulting model can be ex- 
ploited to estimate the wear level. ARMA (Auto Regressive Moving Average) linear models are widely used 



for m odelling stationary time series in general ( Brockwell and Davisl . Il991t Priestley . 1981 : 



Shumwav and Stoffer . 



20061 ) and are very flexible with respect to real- world applications. ICohen and BersteinI (jl99ll) have success- 
fully discrim inated normal and pathological patients using ARMA analyses on acoustic signals from the respi- 
tory tract. IHo and Xie ( 19981) have used a more generalised ARMA model to forecast failure in mechanical 
systems and lEsquef et al.l (|2003l ) have presented a review of acoustic applications of ARMA models, defining 
in detail a protocol f or the analysis of acoustic spectral features. Some ot her examples may be taken from 
the fields of medicine (lA kav. Welkowitz. Semmlow. and Kostia Il99ll ). fina nce (Ghahramani and Thavaneswaranl 



2006 : Wang and W u. 200 8J), languages (Pawlowskii 119971) a nd engineering (Kannatev-Asibu and Dornfeldl (119821) ■ 
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Karialainen et al.l(|2002l) . ILi and Kareeinl(jl990l ). citeground. lOlafsson and Sigbiornssonl(|2006l ). lNa and Rheel(|2002l )). 
Previous research on machining by means of techniques related to time series analysis was conducted by Profes- 
sor Wu and his t e am, and reported in several pap e rs da ti ng back to the la te 1 980s. Amongst oth er we mention 
Kim et all (fl989[ ) . iFassois et al.l (|l989[ ). lYang et al.l (Il985[ ). lAhn et al.l (|l985h and IHand^ dilii) . 

Although few naturally occurring processes are intrinsically linear, the aim of this section of the analysis is to 
understand if ARMA models are suitable for the purpose of representing acoustic emission signals, and if a linear 
approximation could give us enough information about the dynamical process itself. 

Definition 3.1 (ARMA processes) Let us consider a real valued process Xt, < G N. It is called an ARMA {p, q) 
process (combining an AR. {p) and MA (q) model) if 



At = Co 



iXt- 



nX; 



t-p- 



-et + Oiet^ 



+ dqEt-q 



(2) 



where cq is the intercept, p and q are respectively the number of parameters in the autoregressive and moving average 
part of the process. The p + q parameters <j)i ■ ■ ■ (pp^Oi- ■ ■ 9q must be chosen such that \4>i\ < 1, « G {1, • ' ' tP} ^^^^ 
l^jl <l,j (z {I,--- ,q} to ensure stationarity and inv ertibility of the p roce ss. The process et ca l led in novations is 



taken to be a white noise (see e.g. the classical text o uBox et al\ ^20 04 ) or Wrockwell and Davia 1199 A) 



For each level of wear a nd for each recorded time series w e have conducted a full analysis following the Box- 
Jenkins iterative procedure ( Box. Jenkins, and Reinsell . 12004) ). The resulting best model is generally very simple; 
often it has only three or four parameters. We can analyse the results, for example, for an RMS time series 
taken with a fully worn tool. At the right side of figure |3| we can see the autocorrelation function and the partial 
autocorrelation function of the series. Their shapes (the acf 's decay is exponential and the pacf is zero after lag one) 
and an analysis of the spectral dens ity function suggest a simpl e AR(1). Furthermore if we apply the procedure 
in a completely automated manner ( Brockwell and Davisl . Il99ll par. 9.3), the resulting best model is indeed an 
AR(1) model. 



Xt = Co + (t)Xt-i -I- et 

= 7.5499 -H 0.7632At_i -h et 



(3) 



where cg is the estimated intercept, (j) is the sole estimated autoregressive parameter and et are the residuals (i.e. 
the estimated innovations). All the parameters are statistically significant. 

Definition 3.2 (BIC — Bayesian Information Criterion) We define the Bayesian Information Criterion for 
a model with [p + q) parameters and N observations as 



BIC [p, q)^N log &^ (p, q) + {p + q) log N+ 



{N-p-q)\og(l- 



N 



+ {p + q) log 



fe?) 
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where p is the number of parameters in the autoregressive part of the model, q is the number of parameters in the 
moving average part, cP' (p, q) is the residuals variance ca lculated af t er ha ving fitted an ARMA [p, q) model and a1 
is the sample variance of the observations (for details see \Priestleil u98i ) page 375). 
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Figure 3: BIC matrix and autocorrelation functions. 



To better unde rstand the behaviour of the BIC it is useful to give the following representation also presented in 
PriestlevI (|l98l[ ) 

BIC ip,q) = A^loga^ (p^q) + {p + q)logN + o{p + q) (5) 

The BIC must thus be minimised w.r.t. {p, q) in order to select the model which best explains the observations 
with the minimum of parameters. 

At the left side of figure [3] we can see the levelplot for the BIC matrix. For each combination of order {p, q) of 
the AR and MA parts, maximum likelihood estimates for the parameters are computed and BIC(p, g) is recorded 
on the matrix. The minimum BIC value indicates that the best model is that of eq. ([3]) . 

Following the usual procedure to validate the model, we conduct an analysis of whiteness on the residuals. The 
autocorrelation and partial autocorrelation functions for the estimated residuals are both within the white noise 
confidence band for lags strictly positive, and the spectral density function is uniform. Whiteness tests in both the 
time and frequency domains (Ljung-Box for various lags and cumulative periodogram tests) were conducted with 
positive outcome. This model therefore explains the whole linear dependence between the variables in the process. 

At this point we have a model which adequately describes the stationary part of the AE RMS signal. Furthermore 
we can see that the mean of the residuals' variance decreases with increasing tool wear. However the decrease of the 
residuals' variance (with respect to the wear level) is not due to the better explanation of the data by the model, 
but instead to the decreasing number of bursts in the time series. Therefore ARMA models are more suitable for 
the description of the essentially continuous transformed signal part (the one due to plastic deformation). 



4 Point processes 

Taking a closer look at the residuals (figure|4]), we note that bursts are still present, though with a smaller amplitude. 
Therefore the bursts are not explained by the AR(1) model of the preceding section, and therefore not by any 
ARMA(p, g). This indicates a decrease in the mean burst frequency with increasing wear level. If we wish to 
understand the evolution of the acoustic emissions with respect to the wear level of the tool, we must take into 
account the dynamical and statistical process that generates the bursts. 

The acoustic emission bursts are usually associated with micro-fractures in the work tool. They are, effectively, 
singular events (the exponential decay is due to the transducer response). Furthermore the number of bursts seems 
to c hange with t he we a r level. Point proc e sses a re a tool widely used in modelling inherently point phenomena (see 
e.g. IPetri et al.l (|l994l ). iLowen and TeichI ( 20051 )). In this section we consider the bursts as the outcome of a point 
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Figure 4: A typical residuals time series after adapting an ARMA model. We note that bursts are still present. 



process and try to understand the behaviour of this process as the wear level increases. Figure [5] shows the overall 
point process for new tools and the observed waiting times between bursts. The waiting times between events are 
registered and presented in figure [Sl^b) indexed by their order in the sequence of bursts. 

The identification of the events was performed by placing a threshold at various levels to obtain information about 
bursts of different amplitudes. In particular we chose four levels of amplitude 40, 50, 60 and 70. These thresholds 
are expressed in arbitrary units proportional to the signal's amplitude which depends on the data collection chain. 
For an example see fig. [5] where the placement of a threshold at level 20 identifies two burst events in the RMS- 
transformed time series. The waitin g times process is then calculated for all the thresholds considered and is in 
effect a renewal process (jCoxl . Il962l ). It is possible to verify the lack of correlation (second order dependence) in 
the waiting times process by calculating the spectral density function and the complete and partial autocorrelation 
functions. 

The hypothesis of exponentially distr ibuted waiting t imes (Pois son process) and Pareto (power-law) distributed 
waiting times (fractal point process) as in lThurner et al.l (|l997l ) and lLowen and TeichI (|l993[ ) is addressed conducting 
Kolmogorov-Smirnov goodness of fit tests on the p.d.f. for all thresholds and wear levels. The tests resulted not 
significant, so a nother hypothesis f or the distribution must be proposed. 

As stated in iMalevergne et al.l ( 2005 ). the Weibull distribution (also known as Stretched Exponential) can have 
(for typical parameters' values) both the features of an exponential distribution and a Pareto distribution: 
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where /3 and a are, respectively, the scale and the shape parameter {x > 0,a,l3 > 0). 

When a = 1 a Weibull is an exponential distribution with rate A = -i. This kind of distribution, therefore, allows 
us to describe point processes that have waiting times following an exponential behaviour, but that can exhibit a 
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Figure 5: a) The sequence of burst events as a function of time, b) The waiting time between events indexed by 
their order in the sequence of bursts. 



Tabic 1: Maximum Likelihood Estimates and P- Values of Kolmogorov-Smirnov Tests 
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heavier tail. Table [T] summarises the maximum likelihood estimates for the Weibull parameters a and /3. We have 
calculated the estimates for each threshold and each wear level, finding the estimated mean /i and performed a 
Kolmogorov-Smirnov test to see if the Weibull hypothesis can hold. As can be seen from the P-values in the table, 
the test is not significant for each level and each threshold. Table [T] contains only the estimates for the new and 
totally worn tools. The half worn tool has sufficient data for this analysis only with the lowest threshold; the result 
is consistent with the other wear levels (i.e. jj, — 113.15). 

Figure [6] shows the maximum likelihood Weibull fit for the distribution of waiting times between the bursts (new 
tools, threshold equal to 40). The estimated parameters are a = 0.75 and /3 — 83.1. In the inset we note that, for 
all thresholds considered, the estimated mean of the waiting times increases with the tool wear level. For purposes 
of TCM we could therefore monitor the mean of the estimated distribution to decide whether is necessary to change 
the tool or not. 

5 Conclusions 



In this analysis we have shown an application of time series models to acoustic emission signals from metal cutting 
processes. Before this work few stochastic analyses were performed on this particular type of data. After trans- 
forming them by means of the RMS transform, we applied standard time series statistical techniques to explain 
the underlying process that generates the phenomenon under consideration. In particular initially we used linear 
ARMA models together with the well known Box- Jenkins iterative procedure. We found that these linear models 
with a small number of parameters are suitable for the description of the linear contribution of the background part 
of the signal. The variance of the residuals decreases when the wear level increases but this effect is due mainly to 
the decreasing number of bursts in the time series. Even though, for the purposes of TCM, we could in principle 
monitor that variance, we have obtained better results looking directly at the underlying mechanism generating the 
burst process. 

A renewal point process with Weibull distributed waiting times seems to represent the burst process very well. 
In particular, for each threshold and each wear level the Kolmogorov-Smirnov test results not significant. When the 
tool wear level increases the tail of the distribution becomes heavier and the estimated mean of the waiting times 
process increases. The number of bursts is actually correlated with increasing tool wear, therefore we could identify 
the estimated mean as the parameter that should be used to monitor the condition of the working tool (the tool 
will be substituted when the estimated mean becomes sufiiciently high). 

In conclusion, it is important to underline that, from a methodological point of view, both ARMA models and 
Weibull point processes are suitable for modelling the phenomenon (even though ARMA models are somewhat 
limited to explaining the linear contribution of the plastic deformation) . What appears of importance is that the 
shape of the distribution of waiting times changes in the sense that the tail becomes heavier when the wear level 
increases. Therefore the process creating the bursts seems to evolve in a fundamental manner. 

In future work, it would certainly be very interesting to follow the evolution of the distribution shape along the 
whole life of the cutting tool to try to understand better the dynamical properties of the underlying process. 
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Figure 6: Weibull fit for the pdf of the waiting times process. The mean inter-burst waiting time consistently 
increases with increasing tool wear. 
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